Pump-probe spectroscopy of two-body correlations in ultracold gases 
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We suggest pump-probe spectroscopy to study pair correlations that determine the many-body 
dynamics in weakly interacting, dilute ultracold gases. A suitably chosen, short laser pulse depletes 
the pair density locally, creating a 'hole' in the electronic ground state. The dynamics of this non- 
stationary pair density is monitored by a time-delayed probe pulse. The resulting transient signal 
allows to spectrally decompose the 'hole' and to map out the pair correlation function. 

PACS numbers: 03.75.Kk,32.80.Qk,82.53.Kp 



Introduction Bose-Einstein condensation in dilute 
gases is determined by the nature of the two-body in- 
teractions between the atoms [if. These microscopic in- 
teractions manifest themselves in two-body correlations 
and dictate the mesoscopic and macroscopic properties 
of the condensate. Formally, the dynamics of an ultra- 
cold gas is described in terms of field equations For 
dilute gases where only two-body interactions are promi- 
nent, the equation of motion for the field operator that 
annhilatcs (or creates) a particle at position x reads 
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Expectation values of the many-body system can be ex- 
pressed in terms of normal-ordered correlation functions. 
To lowest order, these are the condensate or mean field 
wavefunction, ^(x;t) = (ip(x;t)), the one-body density 

matrix R(x,y;t) = {ip\x;t)ip(y;t)) and the pair corre- 
lation function, $(x,y;t) = (^(x; i)V>(y; t)). For practi- 
cal calculations, the infinite set of equations of motion 
for the many-body problem needs to be truncated. This 
can achieved by expanding the correlation functions into 
cumulants Q. A separation of time or length scales, i.e. 
small collision time vs long free propagation time or small 
effective range of the interaction potential vs large inter- 
atomic distance is required to justify truncation. For 
dilute Bose gases in a macroscopic trap, such an assump- 
tion can typically be made, and within the first-order 
cumulant expansion, the dynamics of pair correlations is 
decoupled from higher order terms Q ■ Alternatively, one 
can work with the correlation functions directly [3j|. In 
both cases, the dynamics of the macroscopic pair correla- 
tion function is described by a Schrodingcr-like equation 
where the mean field enters as a source term |2| or acts 
as an additional potential If we restrict our con- 
siderations to timescales that are much shorter than the 
timescale of the mean field dynamics, the pair correlation 
dynamics are described by a standard Schrodinger equa- 
tion where the presence of the condensate only modifies 



the boundary conditions. The macroscopic pair correla- 
tion function is then given by the two-body wavefunction 
of an isolated pair of atoms, $(r) with r = \x — y\. We 
can thus study the many-body pair correlation dynamics 
by solving the time-dependent Schrodinger equation for 
two colliding atoms 

In the most simplified approach, the effect of the many- 
body pair correlations is captured in a single parameter, 
the scattering length [l| . Measuring the scattering length 
corresponds to an indirect assessment of the pair corre- 
lations. If the two-body interaction is probed in a time 
much shorter than its characteristic timescale, a more 
comprehensive study becomes possible. Here we suggest 
to employ pump-probe spectroscopy to unravel the dy- 
namics of pair correlations in an ultracold Bose gas. This 
requires a combination of ultrafast and ultracold physics, 
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FIG. 1: (Color online) Pump-probe spectroscopy of dynami- 
cal pair correlations: A pump pulse excites population from 
the electronic ground, leaving the pair correlation function 
in a non-stationary state, the 'hole'. A time-delayed probe 
pulse monitors the dynamics of the 'hole'. The orange peaks 
indicate the action where the probe pulse measures pair am- 
plitude. 
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the basic basic feasibility of which has been demonstrated 
in recent experiments on femtosecond photoassociation of 
ultracold rubidium atoms [3, @] • 

Pump-probe spectroscopy Our scheme involves three 
short laser pulses and is sketched in Fig.UJfor the example 
of an ultracold gas of 87 Rb atoms. The pump pulse ex- 
cites population from the electronic ground to an excited 
state, leaving a 'hole' in the initial pair correlation func- 
tion. The 'hole' represents a non-stationary state that 
moves under the influence of the ground state potential, 
cf. Fig. [1] The pump pulse thus induces the dynamics 
of the pair correlations. A time-delayed probe monitors 
these dynamics by measuring the amount of probability 
amplitude in a range of internuclear distances. The mea- 
surement consists of applying simultaneously a photoas- 
sociation and an ionization pulse (combination of red and 
green arrows in Fig. [1]) . The pair density on the ground 
state is thus photoassociated and immediately ionized for 
detection. The largest probe signal is obtained when the 
probe pulse is identical to the pump pulse. The dynam- 
ics are then monitored at the position where the 'hole' 
was created. The spatial region where the probe pulse 
detects pair density is indicated in orange in Fig. [TJ 

For alkali atoms, the initial pair correlation function 
consists of a superposition of singlet and triplet com- 
ponents. The corresponding interaction potentials are 
shown in solid (singlet) and dashed (triplet) lines. For 
clarity's sake, only triplet pair wave functions are de- 
picted in Fig. [TJ Since the respective excited state poten- 
tials differ, the probe pulse is resonant at different inter- 
atomic distances for singlet and triplet. This is indicated 
by the orange peaks representing the probe windows in 
Fig. [JJ (striped for singlet, dotted for triplet). 

Detection of the pair correlation dynamics proceeds 
via the creation of molecular ions and is inspired by 
Refs. 0, ij. In those experiments, the pump pulse 
photoassociates ultracold atoms and the ensuing excited 
state molecular dynamics are detected by the ionization 
pulse. In contrast to that, it is the ground state dynam- 
ics of the many-body pair correlations that are probed 
in our proposal; photoassociation just serves as a means 
for detection. Hence, while the experimental setups envi- 
sioned in this work and realized in Refs. 0,0] are similar, 
the probed physics is rather different. 

Modelling the two-body interaction We consider two 
colliding 87 Rb atoms. Hyperfine interaction couples the 
ground state singlet and lowest triplet scattering chan- 
nels. However, this interaction cannot be resolved on the 
timescales considered below. We therefore assume a su- 
perposition of singlet and triplet components, but neglect 
the effect of hyperfine interaction on binding energies and 
dynamics. The two-body Hamiltonian is represented on 
a grid large enough to faithfully represent the scattering 
atoms. The interaction of the atom pair with the pump 
pulse is treated within the dipole and rotating wave ap- 
proximations. Excitation is considered exemplarily into 



the 0+(5s + 5p 3 / 2 ) and 0~(5s + 5p 3 / 2 ) excited states. The 
pulses are taken to be transform-limited Gaussian pulses 
with a full-width at half-maximum (FWHM) of 10 ps. 
This corresponds to a spectral bandwidth of roughly 
1.5cm" 1 or 45 GHz. Details on the potentials and the 
employed methods are found in Ref. [g]. For a Bose- 
Einstein condensate, a single low energy scattering state 
needs to be considered. The collision energy of this initial 
state is chosen to correspond to 20 [iK with 75% (25%) 
triplet (singlet) character. At higher temperatures, the 
bosonic nature of the atoms can be neglected; and the 
ultracold thermal ensemble is described by a Boltzmann 
average over all thermally populated two-body scattering 
states [7|- 

Modelling the absorption of the probe pulse The dy- 
namics of the non-stationary 'hole' is monitored by a 
combination of a probe pulse and an ionization pulse, 
cf. Fig. [1] This two-color scheme converts the absorp- 
tion of the probe pulse into detection of molecular ions. 
Assuming the probe pulse to be weak and the ionization 
step to be saturated, absorption of the probe pulse can be 
modelled within first order perturbation theory [l(| ■ 
The transient absorption signal is then represented by the 
time-dependent expectation value of a window operator, 

W(r)=^(r p ^ ) 2 e - 2A « 2 ^.A 2 , (2) 

where t p and E P: o denote duration (FWHM) and peak 
amplitude of the probe pulse. p, p is the transition dipole 
moment between ground and first excited state. The 
central frequency of the probe pulse, uj p , determines the 
difference potential, A(r) = V e (r) — V g (f) — hw p . A 
position measurement becomes possible if the difference 
between the ground and excited state potential is suffi- 
ciently large, and, moreover, the spectral bandwidth suf- 
ficiently small to probe only non-zero A(f). Since the 
difference potential vanishes for r — > oo, this implies suffi- 
ciently detuned, narrow-band probe pulses. Fig. [JJ shows 
triplet and singlet window operators (orange peaks) as- 
suming identical parameters for pump and probe pulses. 

Characterization of the 'hole ' As sketched in Fig. [TJ 
the pump pulse carves a 'hole' into the ground state pair 
correlation function. The resulting non-stationary wave 
packet is a superposition of a few weakly bound vibra- 
tional wave functions and many scattering states Q . The 
detuning of the pump pulse from the atomic resonance 
frequency, Al = u>l —ui a t, determines the position where 
the 'hole' is created: For larger detuning, excitation oc- 
curs at shorter distance and populates deeper bound lev- 
els. A pulse energy of £p = 1.5 nJ is sufficient to deplete 
the population within the resonance window of the pump 
pulse (cf. the blue wave function in Fig. [TJ). At higher 
pulse energies, Rabi cycling within the resonance window 
sets in. This leads to more population in the bound levels 
but also to stronger redistribution among the scattering 
states. 
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FIG. 2: (Color online) Probing the two-body correlation dy- 
namics: Absorption of the probe pulse as a function of the 
delay between pump and probe pulses. The pump pulse de- 
tunings are Al == —4.0cm" 1 (a) and Al = —14cm" 1 (b), 
and three different pump pulse energies are shown. 
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FIG. 3: (Color online) Spectra of the transient absorption 
signals shown in Fig. The eigenergies of the interaction 
potential are recovered: Black arrows indicate the position 
of the binding energies of the triplet (upper row) and singlet 
(lower row) levels, grey arrows half-multiples of the binding 
energies. The spectra for 23.5 nJ (blue solid lines) are scaled 
down by a factor of two. 



Dynamics of the 'hole' Fig.[T]also illustrates the time 
evolution of the pair correlation amplitude after a weak 
pump pulse (£p — 1.5 nJ) has been applied: The blue 
curve depicts the wave function just after the pump pulse, 
at t = 24 ps (taking t — to be the time of the pump 
pulse maximum). A probe measurement at that time 
will find no amplitude within the probe window. Due 
to the attractive interaction potential, the 'hole' moves 
toward shorter distances, cf. the brown wave function 
(t = 126 ps). This brings amplitude that is initially at 
larger r into the probe window. Eventually the motion of 
the 'hole' will be reflected at the repulsive barrier of the 
potential. The bound part of the wave packet will remain 
at short distance and oscillate, while the scattering part 
will pass through the probe window once not to return. 

Pump-probe spectra The dynamics of the 'hole' is re- 
flected in the transient probe absorption, i.e. the time- 
dependent expectation value of the window operator (red 
dotted curve in Fig. [2] a): A depletion of the signal due to 
the creation of the 'hole' by the pump pulse, referred to 
as 'bleach' in traditional pump-probe spectroscopy, is fol- 
lowed by a recovery that peaks at 550 ps. At later times 
oscillations due to partial recurrence are observed but a 
full recovery does not occur. Rabi cycling induced by 
larger pulse energies may partially (green solid curve) or 
completely (blue dashed curve) refill the 'hole'. For larger 
detuning, cf. Fig. [2Jd, the 'hole' is created at shorter in- 



teratomic distance, 



48 a for A L 



-14 cm vs 



r ~ 76 ao for = —4 cm" 1 . Obviously, the time to 
move to the repulsive barrier and back is then shorter. 
A faster recovery of the bleach is hence observed - at 
t = 110 ps in Fig. [2b. 



The spectrum of the transient absorption signals shown 
in Fig. [2] can be obtained by filter-diagonalization llj, a 
method allowing to accurately extract frequencies from 
just a few oscillation periods. The spectra are shown 
in Fig. [3J The vibrational energies of the two-body in- 
teraction potential are recovered. Fig. [3] shows further- 
more that Rabi cycling during the pump pulse leads to a 
larger bound part in the 'hole', cf. the increase of spectral 
weights with pump pulse energy. While a direct measure- 
ment of the vibrational populations would be difficult to 
implement experimentally, wave packet spectral analysis 
via probe absorption is fairly straightforward. 

Pure state vs thermally averaged dynamics Pump 
probe spectroscopy of the pair correlations can be applied 
to a condensate as well as a thermal ultracold gas. For the 
timescales considered here, this translates into compar- 
ing the dynamics of a pure state to that of an incoherent 
ensemble. Fig.0]shows the transient probe absorption for 
the two cases. Thermal averaging markedly smears out 
the recovery of the bleach at 550 ps. In principle, ther- 
mal averaging introduces two effects - the finite width in 
scattering energies which is too small to be resolved on 
a nanosecond timescale, and the contribution of higher 
partial waves. The latter becomes particularly prominent 
in the presence of resonances. Fig. 0] compares calcula- 
tions including both singlet and triplet channels (a) to 
those for the triplet component only (b) in order to high- 
light the role of shape resonances. For 87 Rb, shape reso- 
nances are observed at ~160/iK (J = 2) and at ^430 //K 
(J = 6) in the singlet, and for J = 2 at ^290 [iK in the 
triplet channel. The second peak observed in Fig. QJ, at 
1100 ps corresponds to the singlet recovery of the bleach. 
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FIG. 4: (Color online) Absorption of the probe pulse as a 
function of the delay between pump and probe pulses com- 
paring pure state dynamics (black solid lines) to those of a 
thermal ensemble (colored broken lines), (a) Calculation in- 
cluding two scattering channels (singlet and triplet): The 
thermal dynamics are dominated by a singlet shape reso- 
nance while the pure dynamics shows features of both singlet 
and triplet dynamics, (b) Calculation for a single channel 
(triplet): The thermal and pure state dynamics are similar 
with the recovery of the bleach smeared out at higher tem- 
peratures. (Al = —4.0cm _1 , £p — 3.8 nj) 



The shape resonances lead to a much larger weight of the 
singlet contribution in the thermal averages than in the 
pure state s-wave calculation. Since probe absorption in 
the singlet channel occurs at larger distances than for the 
triplet channel, cf. Fig. [TJ the recovery of the bleach is 
observed at later times. This observation opens up the 
perspective of analysing pair correlation functions in cou- 
pled channels scattering near a resonance where tuning 
an external field through the resonance will modify the 
respective weight of the channels. We emphasize that this 
novel pair correlation spectroscopy is possible even in the 
presence of, e.g., three-body losses as long as the decay 
occurs on a timescale larger than a few nanoseconds. 

Mapping out the pair correlation function Pump- 
probe spectroscopy, a well established tool in chemical 
physics, allows for retrieving the amplitude and phase of 
a wave function 



12j . In the present context we can re- 



trieve the pair correlation density operator p(r, r', t) with 
p(r,r';t) = $(r; f)$* (/; t) for a BEC. This is based on 
Eq. |(2J|, where probe absorption corresponds to a posi- 
tion measurement with finite resolution: Different cen- 
tral frequencies, top, define the position that is mea- 
sured, and the difference potential, A(f), together with 
the pulse duration, t p , control the resolution. These mea- 
surements resolve the amplitude of the pair correlation, 
$*(r'; t)\ 2 . The phase information is obtained by chirp- 
ing the probe pulse which corresponds to a momentum 



measurement ■ The window operator then defines a 
finite resolution measurement in phase space [l3| . Col- 
lecting the expectation values for a sufficiently large set 
of window operators with different positions/frequencies 
and momenta/chirps corresponds to quantum state to- 
mography of p{r, r'\ t) Q. 

Conclusions Pump-probe spectroscopy unravels di- 
rectly many-body pair correlations in dilute Bose gas. 
Existing experimental setups 0, H| need to be only 
slightly modified to implement our proposal. In particu- 
lar, transform-limited pulses of about 1cm -1 bandwidth 
are required for detection of the probe absorption via 
molecular ions. Spectral features on a scale of less than 
lem -1 can be resolved for pump-probe delays of a few 
nanoseconds. Pump-probe spectroscopy of the pair cor- 
relation dynamics allows to capture transient states of ul- 
tracold gases such as collapsing condensates. Moreover, 
it can be combined with static external field control. For 
example, tuning a magnetic field close to a Feshbach res- 
onance may enhance the pair density at short and inter- 
mediate distances (ljj . The resulting coupled channels 
pair correlation function can be mapped out despite the 
finite lifetime of the resonance. Future work will consider 
shaped pulses. Once picosecond pulse shaping becomes 
available, the full power of coherent control can be em- 
ployed to study pair correlation dynamics. 
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